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ABSTRACT 



Spitzer Space Telescope infrared observations indicate that several transit- 
ing extrasolar giant planets have thermal inversions in their upper atmospheres. 
Above a relative minimum, the temperature appears to increase with altitude. 
Such an inversion probably requires a species at high altitude that absorbs a sig- 
nificant amount of incident optical/UV radiation. Some authors have suggested 
that the strong optical absorbers titanium oxide (TiO) and vanadium oxide (VO) 
could provide the needed additional opacity, but if regions of the atmosphere are 
cold enough for Ti and V to be sequestered into solids they might rain out and be 
severely depleted. With a model of the vertical distribution of a refractory species 
in gaseous and condensed form, we address the question of whether enough TiO 
(or VO) could survive aloft in an irradiated planet's atmosphere to produce a 
thermal inversion. We find that it is unlikely that VO could play a critical role 
in producing thermal inversions. Furthermore, we find that macroscopic mixing 
is essential to the TiO hypothesis; without macroscopic mixing, such a heavy 
species cannot persist in a planet's upper atmosphere. The amount of macro- 
scopic mixing that is required depends on the size of condensed titanium-bearing 
particles that form in regions of an atmosphere that are too cold for gaseous 
TiO to exist. We parameterize the macroscopic mixing with the eddy diffusion 
coefficient Kzz and find, as a function of particle size a, the values that Kzz must 
assume on the highly irradiated planets HD 209458b, HD 149026b, TrES-4, and 
OGLE-TR-56b to loft enough titanium to the upper atmosphere for the TiO hy- 
pothesis to be correct. On these planets, we find that for TiO to be responsible 
for thermal inversions Kzz must be at least a few times 10^ cm^ s~^, even for 
a = 0.1 /im, and increases to nearly 10^^ cm^ s~^ for a = 10 /im. Such large 
values may be problematic for the TiO hypothesis, but are not impossible. 
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Subject headings: astrochemistry - difFusion - planetary systems - radiative 
transfer - turbulence 



1. Introduction 



In the last tw o decades, we have moved from the first discoveries of planets beyo nd 
our solar system (jWolszczan fc Fraill Il992l : iMayor fc Quelozl 1 19951 : iMarcy fc Butlen Il996l ) to 
having the ability to frame and address questions about the structures of distant worlds. The 
remarkable pace of detections has accelerated to the point that there are now more than ~340 
extrasolar planets currently known, of which more than 55 transit their primariesJl] 



Th e transiting planets, first discovered by iHenry et al.l (l2000l ) and ICharbonneau et al. 
( I2OOOI ). constitute a particularly exciting subsample. Knowing that their orbits are edge- 
on breaks the degeneracy between their masses and orbital inclination angles. Furthe r- 



more , their atmospheres raay be probe d during both primary eclipse (ISeager &: Sasseloy 



2000: 



Hubbard et al 



Rauscher et al 



2007 



2001 



Brown 



Burrows et al 



20011) and secondary eclips e (jSudarsky et al.l 12000 . 



2003 



20081 : iFortnev et all 120081 ). 



Early speculative calculations are now supplemented with a wealth of observational data 
that constrain the atmospheres of several transiting extrasolar giant planets (EGPs). The 
strong doublets at ~5900 A and ~7700 A in neutral sodium (Na) and potassium (K) con- 
tribute significantly to the opacity of planetary atmospheres through much of the opti- 
cal range. The sodium feature has been se en in the transit spectra of both HP 209458b 



Charbonneau et allbooilPesert et al.ll2008h and HD 189733b feedfield et aDl2008h . ICharbonneau et al. 



( I2OO2I ) point out that the depth of the observed Na feature is less t han would be expected 



if neutral sodium were present at solar abundance, liro et al.l (120051 ) suggest that if, as ex- 
pected, HD 209458b is in synchronous rotation, then Na might have condensed into Na2S on 
the colder night side, which could in part explain an inferred low abundance of atomic Na. 



Another possible explanation is a high gray haze (ICharbonneau et al.l l2002l : iFortney et al. 
2003h . 



In the simplest picture of an atmosphere's vertical thermal profile, temperature decreases 
with altitude (with decreasing pressure). In certain cases of high stellar irradiation, how- 
ever, a chemical species that is a strong optical or near UV absorber could, if present at 



^See the catalog at |http: / /exoplanet.eu 
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high altitude, le ad to a thermal inve rsion in which the temperature increases above a rel- 
ative minimum ( Hubenv et al. [2003). Such an inv ersion is well-known in the context of 
the Earth's stratosphere (ISherwood &: Desslerl I2OOII ). where the upper-atmosphere heating 
is caused mainly by UV absorption by ozone. Infrared observations by the Spitzer Space 
Telescope have s uggested that several giant planets have suc h thermal inve rsions, includ 



ing H P 209458b flBurrows et al.ll2007bl: iKnutson et al.ll2008ah . HD 149026b flFortnev et al. 



20061 : iBurrows et all 120081 ) . Tr ES-4 flKnutson et al.l l2009l). XO-lb (IMachalek et all 120081 ) 



and, perhaps, v Andromeda b ( jBurrows et al.ll2008l ). It is not known, however, what species 
is responsible for the additional opacity high in these planets' atmospheres. 

Titanium (Ti) and vanadium (V) both play significant roles in the spectra of low-mas s 
stellar objects, such as M-dwarfs ( iKirkpatrick et al.lll999l : iMartfn et al.lll999l : lLoddersll2002l ). 
Their oxides (TiO and VO) in gaseous form are strong optical absorbers that, if present at 
near-solar abundan ce in the upper atmospheres of highly i rradiated planets, could c ause 
thermal inversions (IHubeny et al.ll2003l : iFortney et al.ll2006l . l2008l : iBurrows et al.ll2008l ). It 
might be difficult, though, to maintain significant quantities of gaseous TiO/VO aloft. First, 
TiO and VO are significantly heavier than the primary constituent of EGP atmospheres, 
molecular hydrogen. In the battle of molecular diffusion against gravitational settling, heav- 
ier molecules settle more strongly and are concentrated at deep layers. Macroscopic mixing 
processes, such as turbulent diffusion or large-scale advective motions, are then required to 
maintain a high abundance of a heavy species at altitude. 



Second, and perhaps more significant, Ti and V might rain out of the upper atmo- 
spheres. In chemical equilibrium, at low temperatur es or high pressures, Ti and V are found 



in a variety of condens ates that form solid grains (IBurrows &: SharpI Il999l : iLodderd 12002 



Sharp fc Burrows! 120071) . Radiatiye equilibrium radiat i ve transfer models of close-in EGP 



atmosphere s (ISudarsky et al.ll2003l : iHubeny et al.ll2003l : IBurrows et al.ll2007al : iFortney et al. 



20061 . 12OO8I ) predict that, in at least some of the planets inferred to have inversions (e.g., HD 
209458b), there is a cooler region below the inversion in which the temperature and pres- 
sure would cause Ti and V to be in their condensed, solid phases. Such a region may aptly 
be called a "cold-trap," in analogy with the Earth's equatorial water cold-trap. Above the 
Earth's equator, the atmosphere cools and reaches a relative minimum at the tropopause, 
the beginning of the stratospheric thermal inversion. As a result, even though water could 
exist in gaseous form at altitude, most water condenses and rains out before ever reach- 
ing t he stratosphere, leaving the stratosphere quite dry (lBrewerlll949l : ISherwood fc Dessler 
2OOII ). If a similar rain-out process occurs in the Ti/V cold-traps of giant planets, their 
upper atmospheres could be similarly depleted of TiO and VO. 



Furthermore, there could be other cold-tra ps. Short-period EGPs are generally presume d 
to be tidally locked in synchronous rotation (iGoldreich fc Pealelll966l : ISpiegel et al.l 120071 ). 
The large contrast between the intense stellar forcing at the substellar point and the min- 
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imal heating at the poles and on the night sides of synchronously rotating planets can 



mospheric general circulation models (GCMs) i 
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the night side of a planet by zon al winds could l e ad to t he condensation and setthng out 
of these odixes, as considered by IShowman et al.l (l2008bl ). Though settling might be slow 
on orbital timescales (depending on condensate size), if the night side remains a sink for 
hundreds of millions of years or more it could inevitably lead to significant depletion of a 
planet's upper atmosphere. 



The same types of mixing processes that might loft heavy molecules to greater altitudes 
than could molecular diffusion might also help loft grains or droplets. If turbulent mixing 
on a macroscopic scale is vigorous enough, the condensates of Ti/V in the cold-trap could 
be stirred up into the hot, more rarified upper atmosphere, where they might reform their 
optically important oxides. 



Nevertheless, recent observational eviden ce suggests t h at th e upper atmosphere of HD 
209458b might indeed be deficient in TiO. iDesert et al.l (120081 ) find that, redward of the 
Na D lines, the near-constancy of the planet's transit radius places a limit on the TiO 
abundance of ~10~^ — 10~^ solar, an amount that, as we show in ^ is insufficient to explain 
the inferred thermal inversion. This observation is not proof that TiO is underabundant in 
the upper atmosphere of HD 209458b, because transit spectroscopy probes only atmospheric 
regions near the terminator. It is conceivable that the terminator is depleted of TiO, but 
the substellar point is not. 



To address this issue, this paper introduces a model for the abundances of TiO and VO in 
the upper atmospheres of highly irradiated EGPs. We conclude that VO does not contribute 
significantly to thermal inversions. We use our model to estimate the minimum amount of 
macroscopic mixing that would be necessary to make the TiO hypothesis viable on several 
highly-irradiated planets. We find that for HD 209458b, HD 149026b, TrES-4, OGLE-TR- 
56b, and WASP-12b, significant macroscopic mixing would be necessary, an amount that it is 
not at all clear actually obtains. Our model predicts that, for TiO to be present in the upper 
atmosphere at sufficient quantity to cause thermal inversion, K^z must be at least a few times 
10^ cm^ s~^ for 0.1-/im particles, increasing roughly linearly with condensate particle size. 
HD 209458b has the most severe cold trap of the planets we considered. There, for particles 
from 0.1 /im to 10 /xm in radius, K^z must be nearly ~10^ cm^ s~^ to nearly ~10^^ cm^ s~^, 
respectively. For larger particles, even greater values of Kzz would be required. WASP- 12b 
has no day-side cold trap in our models; there, we find that K^z must be ~2 x 10'' cm^ s^^ 
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The remainder of this paper is structured as follows: In ^ we present the model. We 
describe how we generate temperature-pressure profiles and spectra of irradiated expolanets, 
and we parameterize the upper atmosphere abundance of a condensing species in the presence 
of turbulent and molecular diffusion. Section [3] contains the results of our calculations, and 
§1] discusses a few remaining complications. Finally, in ^ we summarize our findings. 



2. Modeling the Vertical Distribution of Condensates 



Atmospheres mix their constituents in a variety of ways, at a range of different spatial 
scales. At the micro-level, molecular diffusion operates; at the macro-level, turbulent mixing 
and large-scale advective motions that (including planetary-scale winds) are likely to be 
dominant. How much TiO would have to be mixed through these various processes, up to 
the low pressures of the upper atmosphere of an EGP, in order to cause a thermal inversion? 
To address this question, in § 12. II we produce one-dimensional radiative transfer atmosphere 
models. We then, in § 12.21 examine how titanium condensates are likely to be distributed 
within these atmosphere models and, therefore, how much gaseous TiO is actually likely to 
reach the high parts of the atmosphere, where it would have to be if it indeed is the needed 
extra absorber. 



2.1. The Radiative Model 



To generate temperature-pressure profiles and spectra for the day and night sides of several 



EGPs, we compute model atrn o spheres with th e code CDOLT LUSTY, described inlHubenv et al. 
(l2003h. lSudarskv et all J2003h. burr ows et all fl2006[ ) and burrows et all JpOOsh . This code 
IS an offshoot of the code TLUSTY faubenvl llQSsI : Iflubenv fc Lanzl Il995l ). with chemical 
abundances and molecular opacities appro priate to the (relatively ) cooler environments of 



planetary atniosphe res and brown dwarfs (ISharp &: Burrows! 120071 : iBurrows fc SharpI Il999 



Burrows et al.ll200ll ). Irradiation in planetary atmosphere models is incorporated using Ku- 



rucz model stellar sp ectra, interpolated to temp eratures and surface gravities appropriate to 
exoplanet host stars JKuruczlll97"9l . Il994j . ^003 ). 



In most of our radiative transfer model calculations, we assume gaseous TiO or VO 
throughout the entire atmosphere. We note that analysis based on this assumption is not 
strictly self-consistent, because we use the T-P profiles generated from this assumption to 
calculate where in fact these gaseous species exist. Nonetheless, because our models indicate 
that gaseous TiO and VO have most of their radiative influence high in a planet's atmosphere, 
the model condition that these species are present in cold-trap regions serves mainly to avoid 
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discontinuities in opacity versus depth and, therefore, to aid convergence of the models. This 
procedure should not introduce large errors. 



Burrows et al.l (120061 ) propose a formalism for treating the redistribution of incident stellar 



flux in a planet's atmosphere in which the proportion of day-side bolometric stellar flux that 
is transported to, and reradiated from, the night side is P„. This redistribution parameter 
plausibly ranges between 0, corresponding to all stellar flux being instantaneously reradiated, 
and 0.5, corresponding to the night side receivin g exactly as much ste llar energy as the day 
side (as a result of advective heat redistribution). iBurrows et al.l (120081 ) suggest that in many 
situations P„ is likely to vary between 0.1 and 0.4. In this study, we take P„ = 0.3, and 
apply the redistribution between 0.01 and 0.1 bars, for all models. 



First, we note that it is unlikely that VO plays a crucial determining role in an EGP's 
atmospheric structure. Titanium and vanadium are both trace elements, but Ti is about 
~10 times as a bundant as V: Ti's solar a bundance relative to hydrogen is ~ 10"'', whereas 
V's is ~ 10~^ (lAnders fc Grevessd Il989l ). Furthermore, TiO's optical opacity is generally 
greater than that of VO. Figure [1] presents a comparison of temperature-pressure profiles 
for six models of HD 209458b, including models that contain VO, but not TiO, and models 
containing TiO, but not VO. Although VO does heat the upper atmosphere somewhat, even 
10 times solar abundance of VO is insufficient to cause a true thermal inversion. In contrast, 
just 40% solar abundance of TiO causes a modest thermal inversion, while solar abundance 
of TiO causes a significant inversion in which the upper atmosphere (~millibar) is several 
hundred degrees warmer than the isothermal layer deeper in (~ 1-100 bars). 



We, therefore, neglect the contributions of VO, and investigate what errors in spectra 
and temperature-pressure profiles are likely to be generated by the model assumption that 
gaseous TiO exists throughout the atmosphere. Figure [2] portrays spectra (left panel) and 
temperature-pressure profiles (right panel) for six models of HD 2094 58b. Superpose d on the 
spectrum plot are the four data points for this planet measured by iKnutson et al.l (l2008al ) 
using the Spitzer InfraRed Array Camera (IRAC). One of the models in this figure is a base 
case that has no additional upper atmosphere abso rber (blue curv e ); this model does not 
produce a thermal inversion and fails to match the iKnutson et al.l (l2008al ) data. Another 
model has solar abundance of TiO throughou t the atmosphe r e (red curve); this one has 
the largest thermal inversion and matches the IKnutson et al.l (|2008a( ) data the best. The 
remainder of the models in this figure demonstrate that our assumption that the extra 
absorber is present throughout the atmosphere does not cause large errors either in predicted 
emergent spectra or in temperature-pressure profiles. Included on Fig. [2] is a pair of models 
with 20% solar TiO. One model has this species present from the top of the atmosphere 
down to 0.01 bars, and the other has it throughout the atmosphere. Finally, there is a pair 
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of mo dels with a gray absorbei@ of opacity k = 0.2 cm^ g~^, similar to the Kg of Burrows et al. 
(120081 ) ■ Again, one model has this absorber from the top of the atmosphere to 0.01 bars 
and the other has it throughout the atmosphere. The differences in spectra and profiles 
between the models with absorbers only in the upper atmosphere and those with absorbers 
throughout the atmosphere are minor enough that for the purposes of this simple study we 
proceed with models that have the absorber everywhere. 



Finally, we consider several VO-free models of HD 209458b in order to estimate what 
mixing ratio of TiO would be required to sustain an upper atmosphere inversion that pro- 
duces spectra roughly consistent with measured data. Figure [3] presents spectra (left panel) 
and temperature-pressure profiles (right panel) for five models of HD 209458b, one with no 
TiO, and others with 10%, 20%, 50%, and 10 0% solar abund a nce of TiO. Superposed on the 
spectrum plot are the four IRAC points from lKnutson et al.l (l2008al ). The models with 50% 
and 100% solar abundance of TiO have thermal inversions in the upper atmosphere and are 
consistent with the IRAC 1 (~3.6 /um), 2 (~4.5 /xm), and 3 (~5.8 /xm) points. The fact that 
none of these models matches all of the observed data shows that our theoretical understand- 
ing of radiation from exoplanet atmosp heres is incomplete. Nevertheless, one crucial lesson 
from the data of iKnutson et al.l (j2008bl ) is that since the IRAC 2 point is significantly higher 
than the IRAC 1 point, and since t he photosphere of IRAC 2 is at greater altitude than 
that of IRAC 1 (IBurrows et al.ll2008l ). the planet would seem ineluctably to have a thermal 
inversion, other theoretical uncertainties notwithstanding. Fig. [3] suggests that if TiO is the 
extra absorber its mixing ratio in the upper atmosphere ought to be no less than ~50% of 
the corresponding solar ratio. 



2.2. Modeling Mixing Ratio vs. Altitude 



As Fig. [3] shows, any cold-trap regions cannot deplete the atmosphere of Ti too signif- 
icantly without leaving TiO insufficiently abundant to produce the inferred thermal inver- 
sions. How can we estimate how much the cold-trap region depletes the upper atmosphere of 
TiO? In this section, we describe our model for determining the atmospheric profiles of TiO 
and the amount of depletion in a turbulent cold-trap region. To do so, we introduce the tur- 
bulent diffusion coefficient K^z (IColegroye et al.lll965l : iLewis fc Fegleylll984j : INoU et al.lll988 



Drossart et al.lll990l : iRodrieo et al.lll990h . Kzz parameterizes, in a single number, a variety 
of processes (including turbulence and other forms of macroscopic mixing) that act for each 
species to equalize the number density (or partial pressure) at different spatial locations. 



^The absorber is not truly gray, but rather has an opacity that is a top-hat function between 3 x 10 
and 7 X 10^'' Hz. 
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We start by identifying the cold-trap region(s). On the day side, there are in general 
three regions: [i) the hot, isentropic convection zone, which contains gaseous TiO; {ii) the 
cold trap, which contains titanium condensates; and, [in) the hot, thermally inverted upper 
atmosphere, in which TiO is gaseous. Depending on a planet's temperature-pressure profile, 
it might also have no cold-traps. The phase of titanium compounds as a function of altitude 
in a planet's atmosphere is found by comparing the planet's temperature-pressure profile 
with the titanium condensation curve for the corresponding metallicity. TiO is gaseous 
where the atmosphere is hotter than the condensation curve and titanium is in condensed 
form where the atmosphere is cooler than the condensation curve. These condensation curves 
vary with metallicity. At a given pressure, higher metallicity causes the condensation curve 
to be at a higher temperature. Figure H] illustrates, in the context of HD 209458b, how the 
comparison of profiles to condensation curves yields the location of cold-trap regions. Models 
of this planet with 40% solar TiO and with 100% solar Ti O are presented, along with the 



corresponding condensation curves (jSharp fc Burrows! 120071 ). The points of intersection are 



demarcated, and cold-trap regions are indicated by cyan regions of the profile curves. For 
example, on HD 209458b in the two profiles shown the largest cold-trap region extends from 
roughly 3 x 10^ bars to 10~^ bars. 

In order to understand how much Ti reaches the upper atmosphere of an EGP, we start 
by considering the vertical distribution of different molecular species in an atmosphere with 
a variety of chemical constituents. For a particular molecular species i (where i could be 
TiO, H2, or any other atmospheric constituent), the vertical abundance profile ni[z] varies 
with time as dui/dt = —V-Fi (+ sourcesj — sinksj), where ni[z] is the number density 
of species z as a function of altitude z, and Fj is the flux of species i. If we ignore sources 
and sinks, and express the f lux as a comb inati on of molecular and turbulent diffusion with 



gravitational settling, as per lBrewerl (119491 ) and lChamberlain &: Hunteru (119871 ). then we may 



write the rate of change of the vertical distribution of species i as follows: 

driilz] _ ^ in f^'^i ^ f '91n7' ^ rnig\ ^ 



dt dz { \ dz \ dz kT 



dz \ dz kT J 

In this equation, Di is the coefficient of molecular diffusion, is the molecular weight, 
firrip is the product of the mean molecular weight (/i) and the mass of a proton {nip), g is 
the gravitational acceleration (assumed to be constant), k is Boltzmann's constant, T is the 
temperature, and Kzz is the coefficient of turbulent diffusion in the vertical direction. 

The timescale for achieving steady state in a diffusion problem with edd y diffusion coef- 



ficien t is given by the square of a characteristic length divided by Kzz (iGriffith fc Yelk 



19991 ): 

Tss ~ L'/Kzz . (2) 
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A characteristic vertical length scale is the pressure scale height 

kT 



(3) 



which is not more than ~10® cm (1000 km) for EGPs. Kzz is not well constrained. A va- 
riety of estimates of K^z for Jupiter's atraosphere place it in the range 10^ — 10^ crn.^ s~^ 
teezard et allbood : IPedev fc Loddersll 19941 : lAckerman fc Marlevl[200ll : ISen Jaffel eTaPboOTh 
Estimates i r i the context of brown dwarfs place it in the range 10^—10^ cm^ s~^ ( ISaumon et al. 



20061 . l2007l : iGriffith fc Yellelll999l : iHubenv fc Burrowsl 120071 ). If highly irradiated EGPs are 
at least as vigorously mixed as Jupiter and brown dwarfs are thought to be, then we may 
find an upper bound for t^s by setting K^z ^ 10^ cm^ s~^ in eq. ([2]): 



^ 10^^ S ~ 3 Myr < Tpianet age , 



(4) 



where Tpianet age IS the planet's age. This indicates that, in much less than a planet's age 
(hundreds of megayears to a gigayear or more), the atmosphere has reached a steady state 
and drii/dt = 0. We further assume that the mean net fluxes are zero in the steady state. 
Equation ([T]) may, therefore, be rewritten as follows: 

(5) 



dPj rriig 
dz kT ' 



+ Kz 



dz 



kT 



0, 



Here, the nj of eq. ([T]) has been replaced by the partial pressure of species i, Pi, in accordance 
with the ideal gas law. 



The general solution to this equation is 



Pi[z] = Pi[zo]exp 



dc 



'mi/fimp)Di + Kz 



(6) 



Hp[(\{D, + Kzz) 

where zq is an arbitrary reference height. In general, T, /x, Dj, and K^z are functions of z. 
Notice that eq. ([6]) indicates that if vertical mixing is vigorous {Kzz ^ -Dj), all species will 
have the same vertical scale height Hp. Conversely, in the absence of macroscopic mixing 
{Kzz ^ -D,;), gravitational settling of heavier molecules will cause each species to have its 
own pressure scale height Hp. ^ Hp/{mi/ firrip). 



Since the total pressure varies as 

p[z] = p[zo] exp 



dC/Hp[C] 



20 



(7) 



the mixing ratio of species i, Yi = Pi/P, may be expressed as 

{mi/fimp)Di + Kz 



Yi[z] = yi[^o]exp 
= Yi[zo]exp 



dC 



z 



Hp[C]{D, + Kzz) Hp[C]}\ 
dC {{nii/ firrip) - 1} D. 



Hp\Q 



D^ + Kz 
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Equation ([8]) has the same general form as eq. ([6]), which prompts us to define a "scale height 
of the mixing ratio" Hy^ as follows: 

Hy [z] = Hp[z] ^ r/^'^^^' ■ (9) 



We consider the behavior of this model in cases where is large and where it is small 
relative to Di. If Kzz ^ -Dj, then ify. ^ Hp, which indicates that the mixing ratio of species 
i remains essentially constant with height, as is expected with vigorous mixing. If Kzz ^ -Di, 
the scale height of the mixing ratio Hy^ ~ Hp/ {{mi/firrip) — 1}, which indicates that the 
e-folding height for the mixing ratio of each species is roughly inversely proportional to its 
molecular weight, which we take to be 2.3 in all modelsjfl 

Presumably, the deep interior of an EGP is well-mixed by convection, with Yi constant 
with z for all species i. In the atmosphere, a useful metric of the abundance of a species is 
its mixing ratio relative to its mixing ratio at depth: 

/W - 1|, (10) 

where z = is taken to be the radiative-convective boundary. 



2.2.1. Vertical Distribution of Gaseous Species 



In regions where TiO is gaseous, how does its mixing ratio vary with altitude? In ac- 
cordance with iKittell (119691 ) we estimate the diffusivity of TiO through molecular hydrogen 

as 



D 



TiO 



where A is the mean free path and (f ) is the mean velocity. We estimate the mean free path 
as A ~ l/(ncr), where a is the cross-section for collision. Taking as a typical cross-section 
a ~ 7r(l A)^, and as a typical speed (f) ~ (8A;T/3mTio)"'^''^5 we estimate -DtiO ~ 90 cm^ s~^ 
at 1500 K and 1 bar. Similarly, this argument suggests that at kilobar pressures .Dxio is in 
the range of ~10~^ cm^ s~^, whereas at millibar pressures, it is in the range of ~10^ cm^ s~^. 
This simple es t imate predicts comparable molecular diffusion values to those calculated in 
Rodrigo et al.l (jl990l ) in the context of the Earth's thermosphere. 



•^In the limit of zero macroscopic mixing, even helium will separate from molecular hydrogen. In this 
case, Yi should be replaced by the mixing ratio of species i relative to H2, the background pressure Pi should 
be replaced by the pressure of molecular hydrogen, Ph2j ^-iid /iTOp should be replaced by 2mp. 



We now consider several possible cases of relative magnitudes of molecular and turbulent 
diffusion: 



1. If there is no turbulent mixing, or if K^z ^ -Dtio; eq. ([9]) implies that Hy^-^q ~ Hp/27. 
On HD 209458b, Hp ~ 500 - 600 km. As a result, if K-,z < /^Tio, the scale height 
of the mixing ratio Hy^^^Q ~ 20 km. In the ~14 pressure scale heights from a kilobar 
to a millibar, the mixing ratio would drop by ~370 e-foldings. In short, without 
macroscopic mixing, /xio = at the top of the atmosphere. 

2. At the opposite extreme, if Kzz ^ -DTiO; the mixing ratio of even a heavy molecule 
remains essentially constant so long as the species is in the gas phase. If, for example, 
Kzz = 10^(DTio)5 eq. ([9]) implies that -ffyTio ~ S70Hp. In this case, in a vertical 
distance of ~14 pressure scale heights the mixing ratio of TiO would decrease by only 
-4%. 

3. Finally, consider a regime wherein Kzz is somewhat, but not overwhelmingly, larger 
than -Dxio- If this is the case, /xio might be neither ~0 nor ~1, but rather an interesting 
intermediate value. Indeed, since the mean free path varies inversely with number 
density - and, therefore, D varies roughly with ~ - it seems likely that there 
are some regions of an EGP's atmosphere that are in this regime. For example, if 
Kzz = IO-Dtio at some altitude, then Hyrj^^Q ~ OAlHp. This would imply that the 
mixing ratio drops off with altitude more precipitously than pressure does, but not 
dramatically more. If, at another altitude, Kzz = lOODxiO) then Hy.j,.q ~ 3.7Hp. This 
would imply that the mixing ratio decreases with altitude more gradually than pressure 
does. 



In the general case, where TiO is in the gas phase /xio varies with altitude according to the 
following equation: 

1 rfl-_ r {{mi/firrip) -l}Di 

J.,Hp[Q D, + Kzz • ^ ^ 



2.2.2. Vertical Distribution of Condensed Species 
In the cold-trap region, titanium is in condensate form. Details of the ch emistry b y whic h 



TiO condenses into a variety of Ti-bearing compounds can be found in iLoddersI (l2002l ). 
Because of the complicated nature of the condensation chemistry, it is difficult to specify 
exactly which condensates form and in what quantity. Therefore, we will simply use the 
subscript q to denote the condensates of species i. We assume that the condensates will 
form roughly spherical condensates of radius a. Lacking a first-principles theory for the 
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process of condensate formation, we cannot predict the size-distribution or the modal value 
of particle radiiQ Instead, we use a wide range of values of a that spans several orders 
of magnitude: 0.1 /im < a < 100 fim. To model the distribution of these condensates, 
we start with an expression analogous to eq. ([T]), only for the condensed form of species i, 
taking 71^ to be the number density of condensates of species i. If we consider the cloud of 
condensates to be a gas, it is an extremely rarified one. To a good approximation, it should, 
therefore, behave as an ideal gas. The "partial pressure" of condensates of species i is then 
approximately given by the ideal gas law: Pc^ ~ riakT. 

The argument that the atmosphere has achieved steady state applies here as well. There- 
fore, we proceed in the same manner as before, by setting the vertical flux of condensates 
equal to zero and solving for Pc-. The solution is given by an expression identical to eq. ([6]), 
except with rrii and Di replaced by mc- and Dc^. 

The background pressure is again given by eq. ([7j). If A^^^. = rricjmi is the number of 
atoms/molecules of species i per condensates, then, following eq. ([8]), the mixing ratio in a 
cold-trap region may be written 

P 

= l^^[0]exp 

The ratio of the mixing ratio /j of species i to its value at depth may be written, 

m-^y (14) 



(13) 



Equation ( |T3l) may be simplified by making use of the low Reynolds number expres- 



sion for terminal veloci ty, appropriate for an EGP's atmosphere ( lAckerman &: Marleyl 12001 



Showman et al.l l2008bl ). In this regime, the viscous drag force, balanced by gravity, is 



F„ = Gnrjav, where r] is the dynamic viscosity, a the particle size, and v the speed of 
the particle through the viscous medium. This expression implies a terminal velocity of 

9 r] 

We estimate the molecular diffusivity of condensates Dc^ using the Stokes-Einstein relation: 



^/stokes = n— ■ (15) 



kT 

-Dstokes = • (16) 

birrja 



''Nucleation theory (jPraine fc Salpeteij|l9771 : lDraindll981t ICooper et al.ll2003l ) might at some point pro- 
vide a method of estimating particle size, but the community is currently far from an ab initio theory of 
condensate sizes in a planetary atmosphere. 
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Equation f[T6|) is valid in the low-Reynolds number regime only when the mean free path 
of the molecules constituting the fluid is small compared to the size of the falling body. 
When this ratio, called the Knudsen number (A^Kn = A/a; iKnudsen &: WeberlllQllI ) is not 
small, a correction is needed. The analysis of small particles moving through high l y rari - 
fied media has a long history, dating back to early theore t ical work bv ICunninghanj (jlQlOl ) 
and contemporaneous experimental work ( Millikan. 1913 : Millikan 19231 ). The form of the 
Cunningham-Millik an-Davies "slip factor correction" (DaviesI 1945 ) has been updated 



some- 



what over the years (lEl-Fandylll953l : lBaines et al.lll965l ). We adopt the value from lLi &: Wang 
fl2003h : 

(3 = 1 + iVK,(1.256 + 0.4exp[-l.l/iVKn]) • (17) 

In the nonzero Knudsen number regime, the true terminal velocity and diffusivity are in- 
creased over the Stokes values (eqs. [15] and [16]) by a factor of /3: 

(18) 
(19) 



D 



tokes 
/9-Dstokcs • 



Substituting eq. ([T8]) into eq. ([T9]) . and noting that rric^ = (4/3)7r(acJ^Pci, we have 



(20) 



When eq. (l20l) is substituted into the expression for ify^ , the result is 



H*y^\z] = Hp[z] 



(21) 



Equation ([HD implies that D^, = 4.4x10"^ cm^ s-^/?(a/l/im)'^(r/1500 K)(r//2.5xl0-^ g cm"^) 
If Kzz ^ -Dci and 111^ ^ /x^^p, then 



tokes 

This implies that, at the top of a cold-trap of vertical extent z — Zq, 

L, ^../(/?^^/stokes) ' 



(22) 



(23) 



3. Results 



Titanium oxide, if present in a highly irradiated EGP's atmosphere, can have a dramatic 
influence on both the temperature-pressure profile and the emergent spectrum. Figures [2] 
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and [3] show this in the context of HD 209458b. However, there are other exoplanets where 
stellar irradiation, and, hence, the influence of gaseous TiO, could be even stronger. This 
can be seen in the sharp contrasts displayed between the two panels of both Figs. and [Hi 



enrv et al 



200' 



20051 



Figures and [H] both show models of five planets: 1) HD 209458b (jE 
which receives incident flux of ~ 10^ erg cm~^ s~^; 2) HD 149026b (jSato et al. 
which receives ~2xl0^ erg cm~^ s~^; 3 ) TrES-4 (IMandushev et al.l 120071 ). which receives 
^2x10^ erg cm'^ s' ^ 4) OGLE-TR-56 b Judalski et al.ll2002h . which receives ~6 x 10^ erg cm 



and 5) WASP-12b jHebb et allbood ). which receives ~9xl0^ erg cm ^ s ^ . The mean 



■2 s-i; 



molecular weight /i is taken to be 2.3 for all planets, and logiof? is taken to be 3.00, 3.19, 
2.86, 3.27, and 3.04, for HD 209458b, HD 149026b, TrES-4, OGLE-TR-56b, and WASP-12b, 
respectively (where g is in cm s^^). (See Tabled] for a summary of these values.) 



Figure O presents temperature-pressure profiles for models of these planets with no TiO 
(left panel) and with solar abundance of TiO (right panel). Figure El portrays spectra for the 
corresponding models. We note that, as described in §2.H models with TiO contain solar 
abundance of TiO throughout the atmosphere, including in the cold trap; this procedure 
should not produce large errors. In the profile plots (Fig. [5]) contain condensation curves for 
titanium at 0.32 solar abundance (dashed-dotted black curve), solar abundance (solid black 
curve), and 3.2 times solar abundance (dashed black curve). In the left panel, there is no 
cold-trap, because there is no TiO. In the right panel, as in Fig.|H cold-trap regions are found 
where each planet's profile is on the cold side of the solar abundance condensation curve. 
The solar abundance TiO models for HD 209458b, HD 149026b, TrES-4, and OGLE-TR-56b 
all have cold-trap regions. Since the model for WASP- 12b with solar abundance TiO is never 
colder than the corresponding condensation curve, our model predicts that WASP- 12b would 
not have a day-side cold trap for titanium if it is at or near solar abundance. Figure [6] shows 
that in all five model planets, solar abundance TiO produces a large thermal inversion in 
the upper atmosphere (and cools the lower atmosphere). This is reflected in the spectrum 
as an increase in planet-star flux ratio throughout much of the near infrared. 

How much TiO is likely to survive the cold-traps on these planets and to reach the upper 
atmospheres? In Figs. [7| through [TOl we address this question. 

Figures [TlfTOl all show the results of integrating eq. fl23|) for different assumed values of a 
and Kzz (and for the different T-P profiles of the different planet models). The T-P profile 
(which is determined by finding a radiative equilibrium solution to the radiative transfer 
equation in a plane-parallel atmosphere, as described in §2.11 and in the cited references) is 
related to altitude through the scale height relationship of eq. ([3]). Terminal velocities Vf are 
calculated with eqs. (fT5|) and f|T7|l . 



Figure [7] illustrates how much turbulent mixing is required to achieve nonzero concentra- 
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tions of TiO at the top of the atmosphere of each of the five planets under consideration in 
this paper. This figure presents vertical profiles of f^io various combinations of a (top- 
to-bottom: 0.1 fim to 1 yUm to 10 fim) and K^z (left-to-right: 10^ cm^ s~^ to 10^ cm^ s^^ 
to 10^° cm^ s^^)- Since WASP-12b has no day-side cold trap in our models, its curves are 
independent of a. To achieve nonzero concentrations of TiO at microbar pressures requires 
values of K^z that are very high for a stably stratified region such as the radiative part of an 
EGP's atmosphere {Kzz ^ 10^° cm^ s~^, even for WASP-12b, with no day-side cold trap). 
To achieve a nonzero concentration of TiO at millibar pressures, the requirements on Kzz are 
not quite so extreme, but even still Kzz must be ^ 10® cm^ s~^, even for the smallest particle 
sizes, on all but the hottest planets (OGLE-TR-56b and WASP-12b). At Kzz = 10*' cm^ s-\ 
none of the planets has any appreciable amount of TiO in the upper atmosphere, and only 
OGLE-TR-56b and WASP-12b have nonzero TiO above ~1 bar, even for a = 0.1 /im. 

As can be seen in Figs. [T}13] and Fig.O the thermal inversions reach maximum temperatures 
at roughly millibar pressures. This indicates that /tio ought to be nonzero at these pressures. 
But what nonzero value is required? Because, as Figs. [8] and [9] demonstrate, /tio transitions 
from ~0 to ~1 in a fairly narrow range of Kzz, a reasonable estimate is that, to cause a 
thermal inversion, f^io should be ~0.5 at p ~ lO^'^ bars. 

Another way to visualize the analysis of Fig. [7] is to ask how /tio varies with Kzz at a 
particular pressure level. Figure [8] portrays this relationship at a millibar, on the same five 
planets, for condensate particle sizes 0.1 /im (top), 1 /xm (middle), and 10 /im (bottom). The 
amount of turbulent mixing required to achieve a given /tio tends to increase with particle 
size. The magenta curve for WASP-12b, however, is independent of particle size and is, 
therefore, the same in all three panels. This figure shows that /tio is nearly or nearly 1 for 
most values of Kzz, and makes the transition over roughly an order of magnitude change in 
Kzz- For condensates of radius 10 /im, even OGLE-TR-56b requires Kzz ~ 10^ cm^ s~^, while 
the other three planets with day-side cold traps (HD 209458b, HD 149026b, and TrES-4) 
require Kzz between 10^° and 10^^ cm^ s~^. 

Figure [9] displays how the millibar-level value of Jtio varies with a and Kzz, for the four 
planets with day-side cold traps. The color contours indicate the mixing ratio of TiO at 
10^'^ bars relative to the interior mixing ratio of titanium. The green band indicates the 
combinations required to achieve the fiducial value of /xio ~ 0.5. This figure shows that 
planets transition from /tio ~ to /xio ~ 1 in a fairly narrow range of values of a and Kzz- 
In general, smaller particles and more vigorous mixing produce larger values of /xio at a 
millibar. 

One more way to frame this result is to ask what value of Kzz is required to bring 50% of 
the interior mixing ratio of TiO up to a given pressure. Figure dU] shows precisely this. For 
the four planets with day-side cold traps, color contours indicate the amount of turbulent 
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mixing needed to achieve /xio = 0.5, as a function of a and P. At millibar pressures, each of 
these planets requires from ~10'' cm^ s~^ of turbulent mixing (in the case 0.1-/im particles 
in the cold trap of OGLE-TR-56b) to ~10^^ cm^ s~^ (in the case of 100-/im particles in the 
cold trap of HD 209458b). 

Table [1] also presents the value of the eddy diffusion coefficient Kzz that is required to 
maintain /tio = 0.5 above the day-side cold-trap, for various assumed condensate sizes. 0.1- 
/im particles require from K^z ~ 10'' cm^ s~^ (in the case of OGLE-TR-56b) to ~10^ cm^ s~^ 
(in the case of HD 209458b), and these values increase roughly linearly with particle radius. 
l-fim particles require Kzz ~ 10'' to ^10^*^ cm^ s~^; and 10-/im particles require Kzz ~ 10^ 
to ~10" cm^ s-^ 

Since we lack constraints on both a and Kzz, our uncertainty spans many orders of 
magnitude. Particle radii from 0.1 to 30 /xm or more are not implausible, nor are 
values of Kzz from 10^ to 10^ cm^ s^^, possibly even greater. Only for a narrow range of 
a-Kzz space is the upper atmosphere abundance of TiO at all sensitive to a and Kzz- For 
most of parameter space, there is either not nearly enough upper atmosphere TiO to cause 
thermal inversions, or easily enough. 



4. Caveats 

What should we conclude from the results in In this section, we address a few 
additional complications of our analysis. First, in §4.H we consider how a more sophisticated 
analysis would treat the single-zone model, which considers motions only in the vertical 
direction, presented in ^ Then, in §4.21 we qualitatively describe how horizontal winds, 
which effectively couple the day side of a planet to colder parts (including both the night 
side and the polar regions), influence the day-side upper atmosphere abundance of TiO. 



4.1. Caveats for the Single-Zone Model 

First, the models presented in §3] assume either solar abundance of TiO or no TiO, but 
a planet's interior mixing ratio of titanium might be super-solar. If so, then /tio could be 
lower than 0.5 while still maintaining a 50% solar mixing ratio of titanium in the upper 
atmosphere. If a planet's interior titanium abundance were twice solar, for instance, Jtio 
could be 0.25 for a 50% solar mixing ratio of titanium at the top of the cold trap. Nonetheless, 
the band in a-Kzz space in which / changes by a factor of two from, say, 0.5 to 0.25, is fairly 
narrow in comparison to the range of plausible values. Even an order of magnitude change 
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from 0.5 to 0.05 corresponds to a fairly modest change in K^z of a factor of ~3. 

It is also conceivable that we dismiss too q uickly the possibili t y that VO contributes signif- 



icantly to thermal inversions. According to ISharp &: Burrows! (120071 ). vanadium condenses 



at somewhat lower temperature than titanium. All else being equal, its cold-trap would, 
therefore, be smaller. Still, since even 10 times solar abundance of VO does not produce as 
large a thermal inversion as has been inferred from observations of HD 209458b it is unlikely 
that VO could play a key role in producing thermal inversions. 

There are two ways in which our analysis is not self-consistent, both related to the false 
assumption in most of Figs. [T]l5] that TiO is present at constant mixing ratio throughout a 
planet's entire atmosphere. First, we use the temperature-pressure profiles that result from 
this assumption to find where the cold-traps are. Second, we also use condensation curves 
calculated based on this assumption in finding the cold-traps. In reality, as (the non-self- 
consistently generated) Fig. [7| demonstrates, the mixing ratio of TiO decreases with altitude. 
Lower abundance of titanium condenses at lower temperature for a given pressure. As a 
result, the condensation curves of a self-consistent analysis would have shallower slopes that 
reflect the decrease in mixing ratio of TiO at lower pressures. A more sophisticated model 
would self-consistently take into account the progressive depletion of TiO in calculating both 
T-P profiles and condensation curves. 

Finally, one might ask whether our use of day-side average models could mask the presence 
of local conditions near the substellar point that are hot enough not to have cold traps. For 
the four planets that our model predicts have cold traps (HD 209458b, HD 149026b, TrES- 
4, and OGLE-TR-56b), could the intense insolation near the substellar points raise the 
temperature enough that the cold traps don't exist at these local hot spots? As described in 



Burrows et al.l (120081 ) Appendix D, the average conditions correspond to the ring around the 
substellar point with direction cosine /i = 2/3. However, local conditions at the substellar 
point (where irradiation is 50% greater than the average used) are such that on HD 209458b, 
HD 149026b, and TrES-4, the substellar point flux still receives less flux than the day-side 
average flux received by OGLE-TR-56b, a planet that according to our model has a cold- 
trap. This suggests that even local 2D models of these 3 planets would have no part of the 
day side hot enough to avoid having a cold trap. Moreover, parts of the atmosphere with 
yU < 2/3 would have even larger cold traps than predicted by our models. Even if some 
planets, such as perhaps OGLE-TR-56b, do not have T-P profiles corresponding to cold 
trap conditions at their substellar points, planetary winds would prevent gaseous TiO from 
remaining at such hot locations for long. Coupling between parts of a planet with different 
local conditions is the subject of §4.2[ 
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4.2. Other Depletion Regions 

The actual process by which planetary atmospheres achieve their spatially and temporally 
dependent chemical distributions is more complicated than the model presented above. In 
addition to molecular and turbulent diffusion, planetary-scale winds advect gas and conden- 
sates from the day side to the night side and baclj^, equator to pole and back, and across 
different altitudes. The chemical composition of advected Lagrangian parcels of air can 
change in response to the local temperature and pressure conditions that they experience 
in the course of their circulation. The TiO content above the day-side cold-trap, there- 
fore, depends on a variety of timescales: the timescale for zonal circulation. Tare ~ irRp/u, 
where u is the zonal wind speed; the timescale for gravitational settling of condensates, 
■7"settie ~ ^-^^//^^/stokes' timcscalc for chemical transitions - the formation and destruction 
of condensates, Xchem; and the timescale for heating and radiatively cooling. Trad- 

Consider the influence of the night-side cold-trap. The basic formalism for the vertical 
distribution of a species on the night side (in the absence of advection) is the same as that 
on the day side. Because the night side has its own temperature-pressure proflle (that is on 
the condensed side of the condensation curve of Ti for most of the atmosphere), it will have 
its own (advection-free) vertical distribution of species i, here labeled rii^n- Coupling between 
the day and night sides may be described quantitatively by adding a source/sink term to 
eq. ([T]). This term ought to be proportional to the difference between the number density on 
the night and day sides, and inversely proportional to the circulation time. In the context of 
a simple model, suppose that each 2;-level on the day side is carried by circulatory winds to 
one unique level on the night side (and vice versa). Let Z„^{z) be the monotonic, invertible 
function that maps day-side altitude z to the night-side altitude to which it is coupled. The 
revised version of eq. ([T]), then, is the following: 

dni[z] _ d_ / djM ^ / d\nT ^ rriig^ ^ 



dt dz \ \dz \ dz kT 



dz \ dz kT J 

The corresponding equation for the night side looks nearly identical, except 1) a night-to-day 
coupling function, which is the inverse function of the day-to-night function (Z^ = ^n~^), 
takes the place of Z^; and 2) the sign of the l/rdrc term is reversed. 

For the reasons discussed in §2.21 we still may assume drii/dt = 0. Now, a vertical 



iShowman et al.ll2008bl suggest that the vertical wind speeds (30 m s^^) in their simulations might be 
sufficient to prevent significant night-side settling of condensates smaller than 30 fim on each circulation. 
It is possible, however, that night-side condensate settling over a ~ billion year timescale drains the upper 
atmosphere of TiO. 
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integration yields 



f drii f d\nT mig\ \ ( drii f d\nT iimpg\ 

+ — rKn[^n(C)]-^aC])c^C = 0. (25) 
^circ Jo 

We briefly consider the qualitative properties of this equation. An important unknown is 
the nature of the coupling function Z^. This function depends on the trajectories followed 
by Lagrangian fluid elements that start at different altitudes in the course of their circulation 
around a planet. Although neither entropy nor pressure is likely to be strictly conserved, 
are flow patterns best described as isentropic or isobaric? 

Since a highly irradiated EGP's night side lacks the day side's intense external irradiation, 
its static stability is less and its radiative-convective boundary extends up to lower pressures 
than the corresponding boundary on the day side. Furthermore, the cooler night side has a 
smaller pressure-scale height. For both of these reasons, if the circulation is nearly isobaric, 
ni^-a[Zn{,z)] could actually be greater than rij^ni^;] for z in the day-side cold-trap region. The 
night-side cold-trap could, therefore, actually serve as a source for the day side, and might 
make it easier for TiO to reach the upper atmosphere. 

In contrast, if circulation nearly follows isentropes, parcels on the day side will travel to 
locations on the night side that are a greater altitude above the convective zone than they 
had been on the day side (i.e., Zjyz] > z). In this case, the night side is purely a sink. 



Showman et al.l (j2008al ) define the radiative timescale as the e-folding time for temper- 
ature perturbations (of magnitude 5T) of a T-P profile to decay: Trad ~ 6TpCp/{dF/dz), 
where p is mass density, Cp is the specific heat at constant pressure, and F is the net ver- 
tical flux. Note that in radiative equilibrium dF/dz equals zero. Therefore, the radiative 
timescale as defined above makes sense only in the case of a perturbed atmosphere. In 
regimes of T-P space where Trad is large relative to Tdrc, a fluid parcel will not change its 
entropy by much during the course of its circulation. For HD 209458b, equatorial winds of 
^1 km s~^ imp l y a circ ulation timescale of Tare ~ 3 x 10^ s. According to the analysis of 



Showman et al.l (j2008af ). on this planet Trad > Tdrc at pressures deeper than 1 bar. In the 
deep regions of the day-side cold-trap, therefore, the night side is likely to be sink. If winds 
travel at speeds close to the local speed of sound in the deep atmosphere, shocks might make 
the circulation non-isentropic, despite the long radiative timescale, although we note that 
circulation models generally predict wind speeds well below the sou nd speed at pressures of 



1 bar and greater (IMenou fc Rauscherl l2008l : IShowman et al.ll2008bl ) 



The other influence of coupling to the night side is that parcels of air from the upper 
atmosphere of the day side, above the cold-trap region, will be advected to cooler regions 
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on the night side. If condensates form nearly instantaneously (rchem ^ ^circ), then any TiO 
above the day side's cold-trap condenses and begins to settle while on the night side. On 
HD 209458b, the terminal velocity of condensates is ~4 cm s~^/5(a/10 yum)^. Condensates 
of radius 10 fim, therefore, fall Azgettie = ''"circ/^^/stokcs ~ ^^f^ while on the night side. 

This settling process is countered by turbulent diffusion. The characteristic distance that 
the condensates are lofted by turbulent diffusion in a circulation time is A^turb ~ V^drc-^zz- 
If Azturb < Azsettie, the day side's upper atmosphere is steadily depleted of TiO, as it 
condenses on the night side and settles into the cold-trap region. If A2;turb > ^Zsettie, this 
process might not significantly alter the upper atmosphere TiO mixing ratio. 



5. Summary and Conclusions 

As has been published, an additional upper atmosphere absorber in the optical can pro- 
duce the thermal inversions inferred from observations. There is an oft-quoted hypothesis in 
the literature that the strong optical absorbers TiO and VO are responsible for these thermal 
inversions. Here, we have studied the viability of this hypothesis, with a radiative-convective 
radiative-transfer model and a model of particle settling in the presence of turbulent and 
molecular diffusion. We applied these models to five highly irradiated EGPs: HD 209458b, 
HD 149026b, TrES-4, OGLE-TR-56b, and WASP-12b, parameterizing our results (see Ta- 
ble [1]) in terms of sizes of condensed particles in cold-trap regions and the strength of eddy 
diffusion. Our most important findings are the following: 

• It is unlikely that VO plays a role in producing an upper atmosphere thermal inversion. 

• In four of the five planets considered, a TiO cold-trap is likely to exist between the 
hot convection zone and the hot upper atmosphere on the irradiated day sides of the 
planets. The titanium that is present in such cold-traps is likely to be sequestered in 
a variety of condensates that settle much more strongly than does gaseous TiO. The 
only planet that does not have a day-side cold-trap is WASP-12b, which receives at 
least 50% more irradiation than any other known planet. 

• Macroscopic mixing is essential to the TiO hypothesis. Without macroscopic mixing 
processes, such as turbulent diffusion, a heavy molecular species such as TiO will not 
be present in a planet's upper atmosphere. Although WASP- 12b, for instance, has no 
cold-trap in our analysis, it still requires turbulent mixing of ~10^ cm^ s~^ (see Fig. [H]) 
if TiO is to be abundant above a millibar in its upper atmosphere. 

• Planetary- scale winds that couple the day side of a planet to the colder night side are 
likely to make it even more difficult than is indicated by the models in this paper for 
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enough titanium to reach the upper atmosphere for TiO to produce a thermal inversion. 

Finally, we estimate how much macroscopic mixing is required to loft enough condensed 
titanium above the day-side cold-trap for TiO to cause a significant inversion. If titanium 
is sequestered in condensates of radius a, then our model predicts that, for gaseous TiO to 
be present in the upper atmosphere at sufficient quantity to cause thermal inversion, K^^ 
must have the following values on the following planets: 1) on HD 209458b, K^z must be 
^ 6.2xl0^(a/l /im) cm^ s^^; 2) on HD 149026b, K^z must be 2.4xl0^(a/l ^m) cm^ s"^; 
3) on TrES-4, K^z must be 2.7xl0*^(a/l ^m) cm^ s~^ 4) on OGLE-TR-56b, Kzz must be 
^ 1.2x10'^ cm^ s'^ for a = 0.1 /im, ^ 2.1 xlO*" cm^ s"^ for a = 1 /im, and ^ 8. 7x10"^ cm^ s^^ 
for a = 10 iim; and 5) on WASP-12b, K^z must be ^ 1.6x10^ cm^ s~^ The analysis that 
leads to these estimates neglects the effect of the night-side cold-trap, and, therefore, these 
values should be taken as lower limits. Because both K^z and a are currently unknown, it 
remains to be seen whether TiO can indeed be responsible for thermal inversions in highly 
irradiated EGPs. Though our results suggest that the TiO hypothesis might be problematic, 
they provide a framework in which to assess it, given improved estimates of Kzz and a in 
the future. 
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Table 1. Planetary gravitational acceleration, stellar irradiation, and required K. 



Planet 


' 9 








required K^z 






cm s"^ 


erg cm^ 




a — 0.1 //m 


1 //m 


10 //m 


HD 209458b 


1000 




1.0 


6.2 X 10^ 


6.2 X 10^ 


6.5 X 10^° 


HD 149026b 


1560 




2.2 


2.4 X 10^ 


2.3 X lO'' 


2.6 X 10^° 


TrES-4 


721 




2.4 


2.7 X 10^ 


2.7 X lO'' 


3.0 X 10^° 


OGLE-TR-56b 


1850 




5.5 


1.2 X 10^ 


2.1 X 10^ 


8.7 X 10^ 


WASP- 12b 


1090 




9.3 


*1.6 X 10^ 


*1.6 X 10^ 


*1.6 X 10^ 



Note. — This table gives planetary g and stellar flux {F^), and values of K^z (in cm^ s"-^) 
required to achieve /xio — 0.5 above the cold trap, for particles sizes of 0.1 /xm to 10 /xm. 
The asterisks in the last row are because WASP- 12b has no day-side cold trap. The required 
value of Kzz-i therefore, is independent of condensate particle size. 
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T-P profiles for models of HD 209458b 
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Fig. 1. — Comparison of the influence of TiO and VO. Temperature-pressure profiles for six 
models of HD 209458b are presented, one with no TiO or VO (blue), three with different 
mixing ratios of VO, but not TiO (green), and two with TiO, but not VO. Even 10 times 
solar abundance of VO has less of an influence on the profile than just 40% solar TiO. This 
suggests that it is unlikely that VO plays an important role in determining whether a highly 
irradiated EGP has a thermal inversion. The dimple at ~10~^ bars in the profile for the 
model wit h no TiO/VO is cau sed by the redistribution to the night side, parameterized by 
Pn = 0.3 JSurrows et al.lboooh . 
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Fig. 2. — Comparison of the effect of model gaseous TiO in the whole atmosphere versus gaseous TiO only in the upper 
atmosphere. This figure presents spectra (left panel) and temperature-press ure profiles (right pan el) for six models of HD 
209458b. Superposed on the spectrum plot are the IRAC data points from iKnutson et al.l (l2008al ). There are two models 
with 20% solar TiO, one of which has this mixing ratio throughout the whole atmosphere and the other with TiO added only 
above 0.01 bars. There are two models w i th an absorber, that is gray between 3 x 10^^ and 7 x 10^^ Hz, whose opacity is 
K = 0.2 cm^ g~^ (the k^. of iBurrows et al.l l2007al ): one has the absorber throughout the whole atmosphere and one has the 
absorber only above 0.01 bars. Finally, there is both a model with no TiO and a model with solar abundance of TiO throughout 
the atmosphere. See § 12.11 for a discussion. 




Fig. 3. — Comparison of different mixing ratios of TiO. This figure shows spectra (left panel) and temperature-pressure profiles 
(right panel) for five models of HD 209458b: with no TiO, and with 10%, 20%, 50%, and 100% of solar abundance of TiO. 
Gaseous TiO is assumed throughout the atmospher e, and the redi s tributi on parameter P„ = 0.3 is used for all models. The 
spectrum plot superposes IRAC data points from iKnutson et al.l ( l2008al ). As the abundance of TiO increases, the upper 
atmosphere temperature increases, since it absorbs a greater fraction of incident stellar energy, while the lower atmosphere 
cools because less of the stellar flux penetrates to depth. The models with 50% and 100% solar abundance TiO have significant 
thermal inversions in their upper atmospheres; the models with less TiO do not. It is most significant that higher mixing ratios 
of TiO cause greater planet-star flux ratios over most of the wavelength range. Models with 50% and 100% solar abundance of 
TiO are decent matches to the IRAC 1 (~3.6 /xm), IRAC 2 (~4.5 /im), and IRAC 3 (~5.8 ^m) points, though they do fail to 
match the IRAC 4 (~8.0 /xm) point. However, the 0%, 10%, and 20% solar TiO models, which lack thermal inversions, entirely 
fail to match the IRAC data. 
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Fig. 4. — The TiO cold-trap. Condensation curves (blue) and temperature-presure profiles of 
models of HD 209458b (cyan and magenta) are shown for two mixing ratios of TiO - 40% solar 
(dashed) and 100% solar (solid). The profiles are calculated assuming TiO is present in the 
entire atmosphere at the specified mixing ratio. Magenta (cyan) parts of the profiles indicate 
where the atmosphere is warmer (colder) than the corresponding condensation curve and, 
therefore, where titanium is in the gaseous TiO (condensed) phase. Filled circles (triangles) 
mark the intersections of the 40% (100%) solar mixing ratio profile with the corresponding 
condensation curve. Notice that, at 40% solar titanium, the atmosphere of HD 209458b 
might have two titanium cold-traps, whereas at 100% solar it has just one. 




Fig. 5. — Comparison of temperature-pressure profiles for models of five planets with no TiO (left panel) and with solar 
abundance of TiO (right panel). This figure shows T-P profiles for models without (left panel) and with (right panel) the 
strong optical absorber TiO, for HD209458b (blue), HD 149026b (green), TrES-4 (red), OGLE-TR-56b (cyan), and WASP-12b 
(magenta). The redistribution parameter P„ is set equal to = 0.3. Condensation curves (black) for titanium are superposed, 
showing the locations of the condensation curves at 0.32 solar abundance (dashed-dotted line), solar abundance (solid line), 
and 3.2 times solar abundance (dashed line). The addition of TiO heats the upper atmosphere and cools the lower atmosphere, 
because more of the incident stellar flux is absorbed high in the atmosphere. 
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Fig. 6. — Comparison of spectra for models of five planets with no TiO (left panel) and with solar abundance of TiO (right 
panel). This figure presents the spectra that correspond to the T-P profiles in Fig. [5l for the same five planets. Here, adding 
TiO increases the planet-star flux ratios in most of the IRAC range. In particular, the flux ratios for IRAC 2 (~4.5 fim) and 
IRAC 3 (~5.8 fim) are significantly increased. 
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Fig. 7. — Vertical profiles of /xio for five planets, /xio (abscissa) is shown as a function of 
pressure (ordinate) for different combinations of particle size a and turbulent diffusion coeffi- 
cient Kzz- Left-to-right, Kzz takes on the values 10^ cm^ s~^, 10^ cm^ s~^, and 10"^° cm^ s~^. 
Top-to-bottom, a varies from 0.1 iim to 1 /xm to 10 /xm. HD 209458b, HD 149026b, TrES- 
4, OGLE-TR-56b, and WASP-12b are represented, respectively, with dark blue, green, red, 
cyan, and magenta (same color scheme as in Figs. [5] and [6]). Since WASP-12b has no day-side 
cold trap in our models, its curves are independent of a. At Kzz = 10^ cm^ s~^, none of the 
model planets has any appreciable amount of TiO in the upper atmosphere. 
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Fig. 8. — Kzz required to achieve /xio at p = 10^^ bars, for five planets, for various con- 
densate particle sizes. For condensates ranging from 0.1 /xm (top panel) to 10 /xm (bottom 
panel), /xio is plotted as a function of K^z on each of the five planets considered in this paper 
(same color scheme as in previous figures). The magenta curve for WASP-12b is independent 
of particle size and so is identical in all three panels. A horizontal black line is shown in each 
plot at /xio = 0.5 to aid the eye in identifying the value of K^z that is required to achieve 
this fiducial relative mixing ratio. 
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Fig. 9. — Color-map of /tio at a millibar, as a function of a and K^z on the four planets with day-side cold traps (HD 209458b: 
top-left; HD 149026b: top-right; TrES-4: bottom-left; OGLE-TR-56b: bottom-right). The color contours indicate the mixing 
ratio of TiO at 10~^ bars relative to the interior mixing ratio of titanium, for various combinations of particle size (0.1-100 fim) 
and turbulent diffusion coefficient (10^-10^^ cm^ s~^). The green band indicates the combinations required to achieve the fiducial 
value of /xio ~ 0.5. 
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Fig. 10. — Color-map of log[i^22] (in cm^ s"'^) required to achieve /xio = 0.5, as a function of a and P on the four planets with 
day-side cold traps (HD 209458b: top-left; HD 149026b: top-right; TrES-4: bottom-left; OGLE-TR-56b: bottom-right). 



